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We have calculated the isotropic Cq coefficients characterizing the long-range van der Waals 
interaction between two identical heteronuclear alkali-metal diatomic molecules in the same arbitrary 
vibrational level of their ground electronic state X 1 T, + . We consider the ten species made up of 
7 Li, 23 Na, 39 K, 87 Rb and 133 Cs. Following our previous work [M. Lepers et. al, Phys. Rev. A 88, 

032709 (2013)] we use the sum-over-state formula inherent to the second-order perturbation theory, 
composed of the contributions from the transitions within the ground state levels, from the transition 
between ground-state and excited state levels, and from a crossed term. These calculations involve a 
combination of experimental and quantum-chemical data for potential energy curves and transition 
dipole moments. We also investigate the case where the two molecules are in different vibrational 
levels and we show that the Moelwyn-Hughes approximation is valid provided that it is applied for 
each of the three contributions to the sum-over-state formula. Our results are particularly relevant 
in the context of inelastic and reactive collisions between ultracold bialkali molecules, in deeply 
bound or in Feshbach levels. 


I. INTRODUCTION 

In the field of ultracold matter referring to dilute gases 
with a kinetic energy E = ksT equivalent to a temper¬ 
ature T well below 1 mK, dipolar atomic and molecular 
systems are currently attracting a considerable interest, 
as they offer the possibility to study highly-correlated 
systems, many-body physics, quantum magnetism, with 
an exceptional level of control [1, 2]. Essentially dipo¬ 
lar gases are composed of particles carrying a permanent 
electric and/or magnetic dipole moment, which induces 
anisotropic long-range dipole-dipole interactions between 
particles, that can be designed at will using external elec¬ 
tromagnetic fields. For example dipole-dipole interac¬ 
tions were shown to modify drastically the stereodynam¬ 
ics of bimolecular reactive collisions at ultralow temper¬ 
ature [3, 4]. 

In this perspective the recent production of ultracold 
heteronuclear alkali-metal diatomic molecules in the low¬ 
est rovibronic [5-9] and even hyperfine level [10, 11] is 
very promising. A crucial step common to most experi¬ 
ments is the conversion of ultracold atom pairs into so- 
called Feshbach weakly-bound molecules [12], which are 
then transferred to a desired target state using a coher¬ 
ent laser-assisted process known as Stimulated Raman 
Adiabatic Passage (STIRAP) [13]. This target state is 
often the lowest rovibrational or hyperfine one, but in 
principle it can be any molecular level [14, 15]. Such 
investigations have stimulated a wealth of combined the¬ 
oretical and spectrocopic studies devoted to various polar 
bialkali molecules [16-26]. 

Following Ref. [27] heteronuclear bialkali molecules AB 
are usually identified as “reactive” or “non-reactive”, de¬ 
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pending on whether the reaction 

2AB ->• A 2 + B 2 (1) 

is exothermic or endothermic. This classification, which 
is based on the differences in dissociation energies at the 
equilibrium distance of AB, A 2 and B 2 , is not modified 
by the effect of zero-point energy for molecules in the 
lowest vibrational level [27]. However the energy dif¬ 
ference between the entrance and the exit channels of 
Reaction (1) is so small, that even for molecules iden¬ 
tified as “non-reactive”, Reaction (1) becomes energeti¬ 
cally allowed above a small value ur of the AB vibrational 
quantum number (see Table I). Therefore the initial vi¬ 
brational level can be viewed as a control parameter to 
study ultracold reactive collisions, and in particular their 
statistical aspects [28, 29]. 

In all these collisions, especially in the “universal” 
regime where AB is assumed to be destroyed with unit 
probability, the isotropic Cq coefficient characterizing the 
AB-AB van der Waals interaction —Cq/R 6 (where R is 
the intermolecular distance) turns out to be a crucial pa¬ 
rameter [40, 41]. Therefore, in this paper we compute 
the isotropic Cq coefficient between all possible pairs of 
identical heteronuclear bialkali molecules in the same vi¬ 
brational level of the electronic ground state X 1 E + . This 
represents an extension of our previous article [42], where 
we presented the Cq coefficient between molecules in the 
lowest rovibrational level. We expand the Cq coefficient 
by distinguishing the contributions from transitions in¬ 
side and outside the X state, and we show that the hier¬ 
archy observed in Ref. [42] for v = 0 persists for a wide 
range of vibrational levels. 

Moreover we discuss the validity of the so-called 
Moelwyn-Hughes approximation [43], which aims at cal¬ 
culating in a simple way the Cq coefficient for two 
molecules in two different vibrational levels. We show 
that, to work properly, the approximation must be ap- 
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Molecule AB 

ADo (cm 3 ) 

fR 

A E vr (cm 3 ) 

Refs. 

23 Na I33 Cs 

238 

2 

392 

[30-32] 

23 Na 87 Rb 

48 

1 

212 

[31, 33, 34] 

23 Na 39 K 

76 

1 

246 

[31, 35, 36] 

39 K 133 Cs 

238 

2 

272 

[32, 36, 37] 

87 Rb 133 Cs 

30 

1 

98 

[32, 34, 38] 


TABLE I. Vibrational quantum number dr above which reac¬ 
tion (1) becomes exothermic for the molecules AB which are 
stable against collisions in their v = 0 lowest vibrational level 
(assuming that they are in their lowest rotational level j = 0). 
We compare ADo = 2D tl= o(AB) — D v= o(A 2 ) — D v= o(B 2 ), the 
energy difference between the entrance and the exit channels 
for molecules in the lowest vibrational level, and A E VR = 
2 x [E vr (AB) — _E^ = o(AB)], twice the energy of the heteronu- 
clear molecule in the vibrational level vr. The difference be¬ 
tween these two quantities is the excess energy of reaction (1). 
The references give the ground state potential energy curves 
of AB, A 2 and B 2 , that we used to compute the corresponding 
vibrational spectrum (using the method of Ref. [39]). 


plied separately to each contribution of the Cq coefficient. 
As the approximation requires the isotropic static dipole 
polarizability, we give this quantity for all vibrational lev¬ 
els and all heteronuclear bialkali molecules under study. 
Those calculations are relevant to characterize the colli¬ 
sions between a ground-state and a Feshbach molecule, 
which are expected to limit the efficiency of the STIRAP 
transfer to ground state molecules. 

The paper is outlined as follows: Section II presents 
the theoretical formalism for molecules in the same vibra¬ 
tional level; in Section III we study in details the calcula¬ 
tions for the RbCs molecule, focusing on the convergence 
of our calculations. In section IV we give our numerical 
results in a graphical form for all possible pairs of like 
molecules in the same vibrational level (subsection IV A), 
and we discuss the validity of the Moelwyn-Hughes ap¬ 
proximation (subsection IV B). All those results are pro¬ 
vided as tables in the supplemental material [44]. Finally 
Section V contains conclusions and prospects. The iso¬ 
topes used in this article are 7 Li, 23 Na, 39 K, ®'Rb and 
133 Cs; and we expect isotopic substitution to only in¬ 
troduce minor changes in our results. Atomic units of 
distance (1 a.u. = ao = 0.052917721092 nm)) and en¬ 
ergy (1 a.u. = 2 x Ryd = 219474.63137 cm” 1 ) are used 
throughout this paper, except otherwise stated. Occa¬ 
sionally atomic units will also be used for dipole moments 
(1 a.u.= 2.54158059 D). 


II. THEORY 


is characterized by a coefficient Cq(v) which is equal to 


Ce(v) = 3 


E 

*i V l*2 V 2 


d Y „ ,•/ „/ d\ 


Xi 


4” ^Xv,i 2 v 2 


( 2 ) 


where i[ (i 2 ) denotes the electronic states of molecule 1 
(2) accessible from X through electric-dipole transition, 
and v' x (v 2 ) the corresponding vibrational level. The ma¬ 
trix element of the dipole moment operator char¬ 

acterizes the strength of those transitions for molecule 1, 
and similarly for molecule 2. The prefactor 2/3 expresses 
that the only non-vanishing matrix elements correspond 
to the states | i^v^) and | i' 2 v 2 ) with the rotational quan¬ 
tum numbers j[ = j 2 = 1. The labels ji V 2 V 11 J 2 will 
be most often omitted in the following for simplicity. In 
Eq. (2), A Xv ,i' lV ' x = 1 - E XvO is the energy dif¬ 

ference between the rovibronic levels \i , 1 ,v[ : j , 1 = 1) and 
\X,v\,j\ = 0), and similarly for molecule 2. 

We can separate the contributions from each electronic 
state, namely C e (v) = X/'i' C 1 q 12 (v) where 


crhv)=i-x 


^XVii^Vi ^Xv,i 2 v 2 


IXi; 


■ Al 


( 3 ) 


Following our previous work on« = 0 [42] , we can define 
Cf (u), Cq(v), Cg~ e (v ) and Cg~ e (v) which are related to 
Cf i2 [v) as 

Cl(v) = ct XA=X (v) (4) 

C e 6 (v)= Y, <#’*>) (5) 

i{,i' 2 ^X 

cr g (v) = cr e (v) = y («) > ( 6 ) 

i 2 ^X 

where for Cf _e (i;), we used the symmetry relation 

<7‘ 2 » = Cg ll2 (v). Note that in Ref. [42] addressing 
the interaction between two identical molecules in their 
v = 0, j = 0 level, the definition of Cf -e (v = 0) quan¬ 
tity corresponds to two times the Cg~ e (v = 0) coefficient 
defined in Eq. (6). 

Due to the weak variation of the permanent electric 
dipole moment (PEDM) function with the internuclear 
distance (see Fig. 1 and Ref. [45]) the contributions of 
the vibrational transitions inside the electronic state X 
are very weak, and Cg(v) can be written to a very good 
approximation 


Cl(v) 



( 7 ) 


We consider two identical heteronuclear bialkali 
molecules in the same vibrational level v of the ground 
electronic state A 1 E + . Since we assume that the 
molecules are in the lowest rotational level j\ = j 2 = 0, 
the van der Waals interaction is purely isotropic, and it 


where d v is the PEDM in the vibrational level v related 
to the transition j = 0 —> j' = 1, and B v the rotational 
constant of v, which is almost equal to half of the energy 
difference between the (v, j = 0) and (v,f = 1) levels 
involved in Eq. (3). 
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III. EXAMPLE OF THE RbCs MOLECULE 

In this section we describe in details our calculations, 
by considering the particular case of RbCs, which has 
been widely studied in ultracold conditions both experi¬ 
mentally and theoretically [8, 9, 11, 21, 46-55]. 

One crucial aspect of our calculation consists in col¬ 
lecting potential energy curves (PECs) and permanent 
and transition electric dipole moment (PEDM, TEDM) 
functions. We used up-to-date molecular data extracted 
from spectroscopic measurements available in the litera¬ 
ture, and quantum chemistry calculations from our group 
otherwise (Table II). For instance in RbCs, the PEC of 
the ground state A 1 E + was calculated in Ref. [38] us¬ 
ing the Rydberg-Rees-Klein (RKR) method. The PECs 
for the states A 1 E + , and 6 3 II and their spin-orbit inter¬ 
action are taken from the spectroscopic investigation of 
Ref. [56]. Five additional 1 E + states, i.e. C and (4)-(7), 
and five 1 II states, i.e. B , D and (3)—(5), were also in¬ 
cluded. Their PECs as well as the PEDM of the X state 
and their TEDMs with the X state were calculated in 
our group with a quantum chemistry approach (see Ref. 
[45] for details). If necessary, the experimental or com¬ 
puted PECs are matched to an asymptotic long-range 
expansion following refs [38, 57]. The results for the low¬ 
est electronic states are displayed in Fig. 1. For every 
electronic state, the vibrational continuum is taken into 
account. The energies and wave functions of the bound 
and free levels are calculated using our mapped Fourier 
grid Hamiltonian method [39]. 

For a proper convergence of the calculations, the exci¬ 
tation of the molecular core states must be included in 
Eqs. (5) and (6), beside the one of the molecular valence 
states. As such core states are unknown and quite diffi¬ 
cult to evaluate, we relied on the fact that their excitation 
energy is much larger than the one of the valence state, so 
that we took their contributions into account through the 
excitations of both ionic cores of the diatomic molecule. 
The procedure is described in details in Refs. [58-60] and 
is not repeated here. These contributions are labeled as 
partial coefficients Cf c where the index c refers collec¬ 
tively to core excitations. 

In the rest of this section we will discuss the influence 
of the different electronic states, and of the inclusion of 
the vibrational continua. 

We start by giving our computed values of Cf(u), 
Cg(v) and Cf“» as functions of the vibrational quan¬ 
tum number v (see Fig. 2 and supplemental material 
[44]). The most remarkable trend is the strong decrease 
in the Cg coefficient as a function of v, which follows 
the strong decrease in Cf. This is due to the variation 
of the PEDM, which is maximal in the vicinity of the 
equilibrium distance, and which diminishes at large dis¬ 
tances, and thus for high v. Figure 2 also shows a slight 
enhancement of Cf which, together with the decrease in 
Cf , induces a maximum Cq of 1.494 x 10 5 a.u. for v = 14. 

Near the dissociation limit, the PEDM and corre¬ 
spondingly Cf(i>) tend to zero. The only contribution 


comes from electronically-excited states. The two RbCs 
molecules almost behave like four free atoms, and the Cq 
coefficient can be written 

C 6 {v -A oo) = Cf h ~ Rh + C 6 Cs - Cs + 2Cf h - Gs . (8) 

Equation (8) is actually a very good approximation, since 
our calculations give 22510 a.u., while we obtain 22862 
a.u. by taking the atom-atom coefficient of Ref. [61]. 

Now we examine specifically the contribution Cf of the 
electronically-excited states. On Fig. 3 we present the 

i' i' 

partial coefficients Cg 1 2 as functions of v, for the states 
i[, i' 2 = A, B, C and D. We also include the contribution 
from the excitation of the ionic core states, which we sup¬ 
pose to be independent from v. [Note that the crossed 
contributions i[ ^ i' 2 have to be counted twice to get 

the total Cf (u).] The coefficients Cf 1 * 2 vary significantly 
under the combined effects of the TEDM variations and 
Franck-Condon factors (FCFs). For example, C BB is 
maximum for v = 0 since the equilibrium distances of X 
and B states are very close to each other. Then C BB 
decreases with v due to a significantly worse FCF. The 
latter is still poor in the high-u region, but C BB follows 
the enhancement of the X —> B TDM. On the contrary 

the contributions Cg 1 " 2 coming from the C 1 E + state are 
smaller, since the minimum of X and C states are no¬ 
ticeably shifted, inducing poor FCFs. 

On Figs. 4 and 5 we discuss the convergence of the 
sum-over-state formula (3) with respect to the number 
of included electronically-excited states. Beyond the four 
first excited states ( A , B , C and D) each additional state 
increases the Cf value by about 150 a.u., that is a little 
less than 1%. This gives a total increase of Cf(v) by 
about 1000 a.u. (5%) when we include up to the (7) 1 E + 
and the (5) 1 n states. This increase only represents less 
than 1% of the total Cg for v = 0 (see Fig. 5). To illus¬ 
trate the effects of higher states which are not included in 
our computed values above (see Fig. 2), we can mention 
that the (8) 3 n adds only 0.05% (8 a.u.). Moreover an es¬ 
timate for the other excited states, which have a TEDM 
with X even smaller than (8) 1 n, reveals that these states 
bring a contribution at least 4 times weaker than the state 

( 8 ) 1 n- 

We have also examined the sensitivity of our results 
due to the PECs that we use as input data, in partic¬ 
ular for the A 1 T, + state. As said above, the calcula¬ 
tions of Fig. ?? are performed with a RKR curve, and 
the spin-orbit coupling with the 6 3 n state is also taken 
into account. The Cg coefficients can be evaluated tak¬ 
ing instead the A 1 E + RKR curve without any spin-orbit 
coupling, or a PEC from a quantum chemistry calcula¬ 
tion with or without the short-range repulsive term taken 
from [62]. As shown on Fig. 6, the final results are weakly 
sensitive to the details of the PECs. Even if the differ¬ 
ent PECs support a different total number of vibrational 
levels, the sum-over-state formula (9), which does not fa¬ 
vor any level with respect to its neighbors, washes out 
such differences. The conclusion in this respect would 
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drastically change if we calculate the dynamic dipole po¬ 
larizability for frequencies close to the transition energies 
[60]. 

Now we check the influence of the vibrational continua 
of the electronically-excited states on the partial (Fig. 7) 
and total (Fig. 8) Cq coefficients. For each excited state, 
the influence of the continuum is negligible for the lowest 
v levels, quite abruptly increases for intermediate ones 
reaching up to 50 % of the total C§ value or beyond, and 
drops down for the very last vibrational levels. Since the 
TEDMs are smoothly varying with the distance, these 
dramatic variations are related to the FCFs. For small 
v of the electronic ground state, the wave function is so 
localized around the equilibrium distance, that the FCF 
with excited-states continuum wave functions is vanish¬ 
ingly small. The FCF is enhanced for moderately-excited 
vibrational levels, because a significant part of the cor¬ 
responding wave functions is located around their inner 
turning point, where the overlap with continuum wave 
functions is favorable. On the contrary, for the high¬ 
est v the wave functions are mostly located around their 
outer turning point where the best overlap occurs with 
the bound vibrational levels of the electronically-excited 
states. 

Finally, the last limiting factor for the precision of our 
calculation is the accuracy of the TEDMs which can¬ 
not be easily measured. The comparison of computed 
TEDMs among various methods for the transitions in¬ 
volving the lowest electronic states has been discussed 
for instance in one of our previous work [63], suggesting 
that the observed typical difference of about 3% is repre¬ 
sentative of the TEDMs uncertainty. For the molecules 
where Cf is dominant, the comparison of the calculated 
PEDM for the v = 0 level [44, 45] with reported experi¬ 
mental values is also noteworthy. A recent measurement 
for RbCs [8] reports d v =o = 1.225(3)(8) D, compared 
to the present computed one, d v —o = 1.246 D; a simi¬ 
lar agreement is observed for NaK calculated at 2.783 D, 
compared to an experimental value of 2.72(6) D [64]. It 
is thus difficult to yield a proper value of the uncertainty 
of the computed Cq coefficients, but we conservatively 
estimate it to about 15%. 


IV. RESULTS FOR ALL HETERONUCLEAR 
BIALKALI MOLECULES 

A. Molecules in the same vibrational level 


molecules, we plotted the partial Cf and total Cq coef¬ 
ficients as functions of the binding energy. We see very 
similar trends with respect to the binding energy, hence 
to the vibrational quantum number v. The coefficient Cf 
increases with v except near the dissociation limit (see 
Fig. 9), while the total C§ coefficient strongly decreases 
with v except for LiNa and KRb (see Fig. 10). This is 
due to the same reason than the one invoked in Ref. [42] . 
For LiNa and KRb, the contribution Cf dominates over 
Cf due to the weak PDM; on the contrary for the other 
molecules, although Cf shrinks following the decrease of 
the PEDM with v (see Eq. (7)), it remains the largest 
contribution except for the very last vibrational levels. 


B. Molecules in different vibrational levels 


Now we turn to the calculation of the coefficient 
Cg{v i,v 2 ) describing the interaction between two ground- 
state molecules in two different vibrational levels v\ and 
V 2 (still with j\ = j '2 = 0), either close to or far from 
each other. This is for instance relevant for the inter¬ 
action between a ground-state and a Feshbach molecule. 
The coefficient Cq{v\,v 2 ) reads 


C 6 (vi,v 2 ) = - Y 

v\ iLvL 


d 2 d 2 

AIll Uj Xv 2 ii'2 v 2 




Xv2,i' 2 ^2 


(9) 


where A xv k ,i’ k v' k ^ i k< v kd' k =' i - ^ x > v kdk=o an 6 dxv k ,i k v k 
are respectively the transition energies and transition 
dipole moments for molecules k = 1,2. 

Giving a numerical value Cq(vi,v 2 ) for each couple 
(v\,v 2 ) and for each molecule would be particularly cum¬ 
bersome. Therefore it is appropriate to express the 
coefficient Ce(vi,v 2 ) as a function of the coefficients 
Cg(y 1 ) and Cq{v 2 ) between two identical levels, using the 
Moelwyn-Hughes (MH) approximation [43] 


/-iMH 

°6 


(ui,u 2 ) 


+ 


q(u 2 ) 

a(vi) 

a{vi) 

a(v 2 ) 


C 6 (v 1 ) 
C e (v 2 ) 


( 10 ) 


where a(v) is the static dipole polarizability in the vibra¬ 
tional level v. For polar molecules a(v) can be expanded 
as 


a(v) = a g (v) + a e (v), (11) 


The same method has been applied to all the other het- 
eronuclear bialkali molecules. The input data used for the 
calculations, especially the PECs, are presented in Table 
II. In Ref. [42] we already gave our results for v = 0 
and we compared them with the literature [41, 50, 65- 
67]. But to our knowledge the present article is the first 
study concerning vibrationally-excited levels. The results 
are given on Figs. 9 and 10 and in the supplemental 
material [44]. In order to provide graphs with all the 


where a g (v) and a e (v) respectively denote the contribu¬ 
tions from the transitions inside the ground state X , and 
to the electronically-excited states. In analogy to Cf (u), 
a g (v) is dominated by the strong purely rotational tran¬ 
sition, and thus can be written, to a very good approxi¬ 
mation, 



( 12 ) 
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In order to discuss the validity of the MH approxima¬ 
tion (10), we assume that the contributions from the 
electronically-excited states can be reduced to a single 
effective transition towards (i*,v*), which means 


a e (v) 


9(l 2 

z a Xv,i*v* 

3 ^Xv,i*v* 


(13) 


and 


C e 6 (v) = 


C e R 


= c§~ e 


(v) = 


rl 4 

a Xv,i*v* 

3Axv,i*v* 

3(2 B v + A x v ,i*v) 


(14) 

(15) 


Now we can compare the Cq coefficients between two 
different vibrational levels v\ and v 2 , calculated either 
directly or with the MH approximation. The direct cal¬ 
culation gives for each contribution 


Cq (vi,v 2 ) = C$(vi,V 2 ) + C e 6 (v ly v 2 ) 

+ C 6~ B (v 1,V 2 ) + Cl~ S (V!,V 2 ) 
d 2 d 2 

C§(v 1 ,v 2 ) = 


3(B V1 + B V2 ) 


Cl(v u v 2 ) = 


Cl e (vi,v 2 ) = 


cr e (v x,v 2 ) = 


9d 2 d 2 

,i^v^ Uj Xv2 ,^2 v 2 

H - ^ Xv2,i2 V 2 ) 
‘2‘dv 1 d 2 Xv 2 ,i*vZ 
3(2 B Vl + Ax« 2 ,»*u*) 

2d 2 d 2 

zu Xv 1 ,qv‘ u v 2 

3(A x vi ,i*v* +2 B V2 ) 1 


(16) 

(17) 

(18) 

(19) 

( 20 ) 


C^ n {vi,v 2 ) ft dl 1 dl 2 /3(B Vl + B V2 ). In the latter case, 
the PDM d Vk , as well as a g (vk ) and Cl(vk), tend to 
zero, and so we obtain Cq(vi,v 2 ) ft C^ r (vi,v 2 ) ft 

‘2d’Xv 1 ,ilvld J Xv2,i2vz/3{Ax v i,ilvl + ^Xv 2 ,i 2 v^)- 

Those analytical estimates suggest that it is more ap¬ 
propriate to apply the MH approximation (10) to the 
partial coefficients Cf (rq, ^ 2 ), Cl(vi,v 2 ) Cl~ e (vi,v 2 ) and 
Cl~ s (vi,v 2 ) separately, rather than to the total coeffi¬ 
cient Cq(vi,v 2 ), namely 




tcr^iy 1 






a e {vi) 

“g(«i) 

a-ejv 2 ) 
a e (fi) 

Q! e (ui) 

a e (v 2 ) 
a e (v 2 ) 


C|(vi)] 




M 


Cl(v 2 ) 

Cliv,) 

c e 6 (v 2 ) 

C§(v 1 ) 


a g( v l) 


Ue(.V 2 ) ~ 6 


C e 6 (v 2 ) 


QgM 

d!e(wi) 

Q!e(^l) 

a g {v 2 ) 


Cl(vi) 

Cl(v 2 ) 


( 21 ) 


( 22 ) 


-1 


(23) 


-1 


-1 


(24) 


Note that Eqs. (19) and (20) are not equivalent when 
vi ^ v 2 . In order to use the MH approximation, we 
need to apply Equation (10) using the following ingredi¬ 
ents: a(vk) is obtained by applying Eqs. (11)-(13) with 
v = Vk (k = 1,2); and Cg(vk) is obtained by adding 
Eqs. (7), (14) and (15) for v = Vk ( k = 1,2). It is then 
straightforward to see that the value of C^ m {v\,v 2 ) ob¬ 
tained in this way cannot be equal to Cq (iq, v 2 ) given in 
Eq. (16). However, to go beyond this general statement, 
it is worthwhile to look closely at some particular cases. 

For weakly polar molecules like LiNa or KRb in 
the ground vibrational level, Cf{v = 0) and Cq(v = 
0) are comparable, although the ground-state polar¬ 
izability is dominant, a g (v = 0) 3> a e (v = 0) 
(see respectively Tab. II and Fig. 1 of Ref. [42]). 
If we compare Eqs. (9) and (10), by neglecting a e 
(see Eq. (13)), we see after some calculations that 
the MH approximation is not valid. On the con¬ 
trary we can show that it works properly either for 
strongly polar or vibrationally highly-excited molecules. 
In the former case, which corresponds to all the 
molecules in low-lying vibrational levels except LiNa 
and KRb, the purely rotational transition is dominant, 
i.e. a g (v k ) > a e (vk) and Cl(v k ) > C§(v k ), Cl~ e (v k ) 
for Vk = 1 and 2, and so we obtain Cq{v i,v 2 ) ft 


This is confirmed by the results of Fig. 11, which displays 
the total Cq coefficient between two RbCs molecules, one 
being in the ground rovronic level, and the other in an 
arbitrary vibrational level v. As long as the Cf contribu¬ 
tion is dominant, we see that is sufficient to apply the MH 
approximation to the total Cq coefficient. But when the 
contribution of the electronically-excited states become 
significant, i.e. from v ft 80 (see Fig. 2), it is necessary to 
apply the MH approximation to each contribution sep¬ 
arately. This requires to know the coefficients C|(ufe), 
C§(vk), and the static dipole polarizabilities a g (vk ) and 
a e (vk), for the two levels of interest V\ and v 2 . To that 
end we give in the supplemental material those five quan¬ 
tities for all vibrational levels v and all molecules [44], 

V. CONCLUSIONS AND OUTLOOK 

In this article, we have computed the isotropic Cq(v) 
coefficient between all pairs of identical alkali-metal di¬ 
atomic molecules made of 7 Li, 23 Na, 39 K, 87 Rb and 
133 Cs, lying in the same vibrational level v of the elec¬ 
tronic ground state A 1 E + . Following our previous work 
[42] we have expanded C e (v) as a sum of three con¬ 
tributions, Cq(v), Cg(v) and C| _e (u), coming respec- 





















6 


tively from transitions inside the electronic ground state, 
from transitions to electronically-excited states, and from 
crossed terms. To a very good approximation, Cj?(v) is 
dominated by the purely rotational transition | X,Vk = 
v,jk = 0} —> \X, v k = v,j k = 1) (k = 1, 2), and so follow¬ 
ing the permanent electric dipole moment d v , it strongly 
decreases with the vibrational quantum number v. In 
comparison the variation is Cq(v) is much smoother, and 
CT» is always vanishingly small [44]. 

We observe that the hierarchy between contributions 
that we established in our previous paper for v = 0 is 
observable for a wide range of vibrational quantum num¬ 
bers v. For all molecules but LiNa and KRb, Cf(u) is 
very strong, and so we expect the effect of mutual orien¬ 
tation described in [42] to occur below the intermolecular 
distance R* = / By ) 1 / 3 . In all other cases (LiNa, KRb 
and other molecules close to the dissociation limit), Cf (i>) 
is dominant, and there is no mutual orientation. 

In order to characterize the van der Waals interaction 
between molecules in different vibrational levels v\ and 
V 2 , we also discuss the validity of the Moelwyn-Hughes 
approximation, in which the coefficient Cq(v 1 ,^ 2 ) is ex¬ 
pressed as a function of C 6 (v 1 ) and Cq{v 2 )- We show that 
it is more appropriate to apply the approximation to the 
partial coefficients Cf(vi,w 2 ), Cq{vi,V2), Cg~ e (v i,i> 2 ) 
and Cg~ s (v 1 , 1 ) 2 ), than to the total one C§(v\,V 2 )- Such 
conclusions can be extended to the interaction between 
two (different) molecules in arbitrary rovibrational levels. 


In particular for rotationally-excited levels, anisotropic 
C|, Cr and Cg~ s coefficients would come into play, 
which could also be calculated with the Moelwyn-Hughes 
approximation. Another readily accessible extension con¬ 
cerns the long-range interaction between a molecule in 
their lowest v = 0, j = 0 level, and a Feshbach molecule. 
Indeed we have shown [68] that the dipole polarizability 
of a diatomic molecule in a weakly-bound Feshbach level 
is well approximated by the one of the uppermost level 
of a single potential curve, itself very close to the sum of 
the two individual dipole polarizabilities [44]. 

Appendix A: Molecular data 

In this appendix we specify the input data - potential- 
energy curves, permanent and transition dipole moments 
and spin-orbit couplings - that we used in our calcula¬ 
tions. In Table II we give in particular the references for 
available experimental data. All the quantum-chemical 
data that we used were computed in our group [45, 60]. 
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Molecule 

experimental 

state 

Experimental excited 
l states 

Long 

Range 

ab initio PECs [60] 

SOCME 

PDMs 

TDMs 

Notes 

KCs 

[37] 

A 1 E + , fdn [69], B 1 !! 
[70], E 1 Y.+ [71] 

[69, 71, 72] 

(3,5,6) 1 E + , (2, 3) A n 

(b/A) [69] 

[45] 

[60] 

B l U : RKR + ab initio + long 
range. 

E 1 E + : modified empirical PEC 
to converge toward the 
unperturbed atomic limit (4s+5d) 

KRb 

[73] 

A 1 E + [74], b J II [75], 

C 1 E+ [76], B 1 n [77], 
D 1 n [78], (3) 1 n [76] 

[72] 

(4,5) I E+ 

(P/A) [18] 

[45] 

[60] 

A 1 S' 1 ", 6 d II : spectroscopic data 
on a very limited range =>• shift of 
ab initio PECs to fit with data. 
SOCME from NaRb is used. 

(1,2) 1 II, C ,1 E'*" : RKR + ab initio 
+ long range 

(S) 1 !! : ab initio + experimental 
T e 

RbCs 

[38] 

A 1 E + , b 3 n [56] 

[38, 72] 

(3 - 7) 1 £ + , (1 - 5)^n 

(b/A) [56] 

[45] 

[60] 


LiNa 

[79] 

A X E+ [80], CPS’ 1 - [81], 
B 1 E+ [82] 

[72, 79] 

(5) 1 E+, (1, 2, 3) 1 11 

— 

[45] 

[60] 

Weak SO interaction not 

included in the calculations. 

LiK 

[83] 

A X E+ [84], C 1 E+ [85], 
B 1 !! [85] 

[72, 83] 

(4,5) 1 E+, (2, 3) x n 

(b/A) [18] 

[45] 

[60] 

Rescaled SOCME from NaRb is 
used. 

A 1 E + , c 1 E+, B 1 n : RKR + ab 
initio + long range 

LiRb 

[86] 

7FTT [87], [87], 

D l U [87] 

[72, 86, 87] 

(2,4, 5) 1 E + , (3) x n 

(b/A) [18] 

[45] 

[60] 

SOCME from NaRb is used 

LiCs 

[88] 

B L n [89] 

[72, 88, 89] 

(2 - 5) I E+, (2, 3) x n 

(b/A) [19] 

[45] 

[60] 

SOCME from NaCs is used 

NaK 

[35] 

6 3 n [90], B 1 !! [91], 
C 1 E+ [92], d 1 n, d 3 n 
[93] 

[35, 72, 93] 

(2,4,5) 1 E+, (2 -4) 1 n 

(d/D) [93], 
(b/A) [18] 

[63] 

Rescaled SOCME from NaRb for 
(b/A) is used (qualitative 
agreement 6% with ab initio from 

[941) 

NaRb 

[33] 

a 1 e + , fdn [is], s 1 n 

[95], C 1 E+ [96], D 1 n 
[97], (4) 1 IT [98] 

[18, 33, 72, 
95-97] 

(4 - 5) 1 E+, (3, 5) 3 n 

(b/A) [18] 

[63] 


NaCs 

[30] 

a 1 e + , fdn [ 19 ], s 3 n 

[99], (3) 1 IT [100] 

[30, 72, 99, 
100] 

(3,4,5) 1 E + , (2,4,5) 1 n 

(b/A) [19] 

[63] 



TABLE II. References for all the electronic PECs used in the calculations (SOCME = Spin-Orbit Coupling Matrix Elements). 
For heavy molecules (RbCs, KRb, KCs) we have also added the short range core-core repulsion from [62] to the ab initio 
PECs. This term is less important for lighter molecules. For the ab initio PECs of LiRb and NaRb we have extrapolated the 
short-range potential by comparing our ground state ab initio PEC to the RKR potential, for the other molecules the term is 
not included. 



FIG. 1. (Color online) Selected permanent and transition dipole moments (a) and potential-energy curves (b) for the lowest 
singlet states of RbCs: A 1 E + (solid line), A 1 E + (dashed line), B 1 !! (dash-dash-dotted line), C 1 E + (dash-dot-dotted line), 
D 1 II (dash-dotted line). The curve labeled X-X corresponds to the permanent dipole moment of the A' 1 E + state, and the 
others to transition dipole moments from the X state. 
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FIG. 2. (Color online) (a) Permanent electric dipole moment d v of the 8 'Rb 133 Cs vibrational levels v in the frame fixed to the 
molecule (see Eq. (7)). (b) Total C& coefficient (solid lines) as well as Cf (dotted lines), Cq (dashed lines) and Cf -e (dasli-dotted 
lines) contributions for the interaction between two RbCs |X 1 E + ,v,j = 0) molecules as functions of the vibrational quantum 
number v. The Ce value indicated by the arrow corresponds to the last vibrational level v = 137. 



FIG. 3. (Color online) Main partial coefficients Cg 1 * 2 ( v ) (see Eq. (3)) as functions of the vibrational quantum number v, related 
to the electronically-excited states A 1 E + (including spin-orbit coupling with b 3 II), B 1 !!, C 1 E + and D 1 !! of RbCs, and to the 
excitation of the core states (see text), labeled with the index “c”. The Ce c (v) and C<p c (v) which are below 200 a.u. for all v 
are not plotted, for clarity. Note that the range of the y axis is split into two parts with a different scale for convenience. 
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FIG. 4. (Color online) Convergence of the C| coefficient between two RbCs \X, v,j = 0) molecules with respect to the number 
of electronically-excited states i' max included in Eq. (3), as a function of the vibrational quantum number v: (i) i' max = 4, (ii) 
ijnax = 6, (iii) = 8, (iv) i' max = 10, (v) i [ nax = 11. The curve (vi) corresponds to a calculation including 11 states without 
including the excitation of the ate”"” 



FIG. 5. (Color online) Same as Figure 4 for the total Ce coefficient in RbCs. All curves for cases (i) to (v) are almost 
indistinguishable at this scale. The coefficient CS(v) is also shown. 



FIG. 6. (Color online) Influence on the coefficient Cg of the origin of the RbCs A 1 E + potential-energy curve, as a function of 
the vibrational quantum number of the X state. The dotted lines correspond to the RKR A 1 E + potential coupled by spin-orbit 
interaction to the RKR 6 3 I1 potential; solid lines correspond to the RKR A 1 E + potential without spin-orbit interaction; dashed 
and dash-dotted lines correspond to ab initio potentials, respectively with and without short-range repulsive wall. 
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FIG. 7. (Color online) Influence on the partial coefficients Cg 1 of the vibrational continua of the electronically-excited states 
i' = A (without symbols), i' = B (circles), i' = C (squares), i' = D (diamonds), as a function of the vibrational quantum 
number of the A' state. Solid lines: vibrational continua included. Dashed lines: vibrational continua not included. 



FIG. 8. (Color online) Same as Figure 7 for the Cg (circles) and Cg (no symbols) coefficients. 



FIG. 9. (Color online) Partial coefficient Cg for the interaction between two identical heteronuclear |A' 1 S + ,n,ji = 0) molecules 
as a function of the binding energy of the vibrational level v. 
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binding energy (cm *) 


FIG. 10. (Color online) (a) Permanent dipole moment d v of the vibrational levels of the ten species of bialkali molecules in 
their own frame (see Eq. (7)). (b) Total Cg coefficient for the interaction between two identical heteronuclear |A' 1 E + ,'U,j = 0) 
molecules as a function of the binding energy of the vibrational level v. 



FIG. 11. (Color online) Total C6 coefficient for the interaction between one RbCs molecule in the ground rovibronic level 
lA' 1 !!" 1 ", vi = 0 ,ji = 0), and another one in an arbitrary vibrational level \X 1 T, + ,V2 = v,j2 = 0), as function of the vibrational 
quantum number v. The solid line corresponds to the numerical calculation; the dashed line corresponds to the Molwyn-Hughes 
approximation (10) applied on the total Ce coefficient; and the circles correspond to the Molwyn-Hughes approximation (10) 
applied on each contribution Cg, Cg, Cf _e and Cg S separately using Eqs. (21)-(24). 





















